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ABSTRACT 

We  examine  the  impact  of  incorrect  atmospheric  correction,  specifically  incorrect  aerosol  model  selection,  on 
retrieval  of  bio-optical  properties  from  satellite  ocean  color  imagery.  Uncertainties  in  retrievals  of  bio-optical  properties 
(such  as  chlorophy  ll,  absorption  and  backscattering  coefficients)  from  satellite  ocean  color  imagery  are  related  to  a 
variety  of  factors,  including  errors  associated  with  sensor  calibration,  atmospheric  correction,  and  the  bio-optical 
inversion  algorithms.  In  many  cases,  selection  of  an  inappropriate  or  erroneous  aerosol  model  during  atmospheric 
correction  can  dominate  the  errors  in  the  satellite  estimation  of  the  normalized  water-leaving  radiances  (nLw),  especially 
over  turbid,  coastal  waters.  These  errors  affect  the  downstream  bio-optical  properties.  Here,  we  focus  on  only  the 
impact  of  incorrect  aerosol  model  selection  on  the  nLw  radiance  estimates,  through  comparisons  between  Moderate- 
Resolution  Imaging  Spectroradiometer  (MODIS)  satellite  data  and  in  situ  measurements  from  AERONET-OC  (Aerosol 
Robotic  NETwork  -  Ocean  Color)  sampling  platforms. 

Keywords:  aerosol  models,  normalized  water-leaving  radiance,  MODIS,  AERONET-OC,  atmospheric  correction, 
ensembles 


1.  INTRODUCTION 

The  main  challenge  in  atmospheric  correction  is  the  estimation  and  the  removal  of  the  path  radiance  from  the 
top-of-atmosphere  (TOA)  radiance  values  recorded  by  the  satellite  sensor  [1].  The  path  radiance  contains  both  Ray  leigh 
and  aerosol  scattering  components  and  can  contribute  about  90%  of  the  TOA  radiance  [2,3].  In  the  current  version  of  the 
NASA  ocean  color  atmospheric  correction  processing  code,  there  are  80  aerosol  models  to  choose  from  to  calculate 
spectral  aerosol  radiance  [4].  This  is  a  more  complex  modeling  system  than  Gordon  and  Wang's  previous  12  model  set 
[5].  Based  on  the  spectral  slope  of  the  aerosol  reflectance  in  the  NIR  bands,  the  two  most  appropriate  aerosol  models 
(from  the  entire  set  of  80  models)  are  retrieved  and  used  for  estimation  of  the  aerosol  radiance  in  the  visible 
wavelengths.  Are  we  appropriately  selecting  these  models? 

We  have  tested  all  80  aerosol  models  individually  with  data  sets  collected  from  three  AERONET-OC  sites 
(Venice,  Martha's  Vineyard,  and  Gulf  of  Mexico)  [6,7,8].  First,  we  derive  nLw  from  MODIS  1km  resolution  imagery'  at 
the  locations  of  the  AERONET-OC  sites,  using  the  aerosol  models  selected  automatically  from  the  standard  atmospheric 
correction  scheme  [9].  We  compare  these  satellite  values  to  the  in  situ  measurements.  We  then  reprocess  the  satellite 
imagery'  using  all  80  aerosol  models  individually  and  again  compare  the  products  to  the  in  situ  measurements.  This  is 
performed  in  order  to  determine  the  “optimal”  aerosol  model  for  each  individual  point  at  the  AERONET-OC  location  for 
each  individual  scene  (MODIS  image).  The  optimal  model  is  the  aerosol  model  that  yields  nLw  closest  to  the  in  situ 
values.  We  determine  the  optimal  aerosol  model  at  a  single  wavelength  and  can  either  use  that  model  for  the  remaining 
visible  wavelengths  or  determine  new  optimal  aerosol  models  at  each  remaining  visible  wavelength. 

Section  2  (Atmospheric  Correction  Using  Aerosol  Models)  briefly  describes  how  aerosol  models  are  used 
during  atmospheric  correction.  This  section  also  describes  the  set  of  aerosol  models  used,  how  those  models  are  chosen 
during  atmospheric  correction,  and  possible  errors  associated  with  the  models.  Section  3  (Automatic  vs  Optimal  Model 
Selection)  compares  the  result  of  standard,  automatic  selection  of  aerosol  models  to  the  optimal  model  selection.  This 
section  shows  that  out  of  the  possible  80  aerosol  models  to  choose,  there  is  most  likely  a  set  of  models  that  w  ill  produce 
nLw  values  that  closely  match  the  in  situ  values.  However,  during  standard  atmospheric  correction,  we  are  usually  not 
properly  selecting  these  models.  This  paper  concludes  with  Section  4  (Conclusions)  and  Section  5  (Future  Work). 
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2.  ATMOSPHERIC  CORRECTION  USING  AEROSOL  MODELS 


The  top  of  the  atmosphere  (TOA)  radiance  is  defined  as  follows: 

Lt(X)  =  LXX)  +  La(X)  +  Lra(X)  +  t(X)  LwC(X)  +  T(X)  Lg(X)  +  t(X)  to(X)  cos0o[Uv(X)  ]N,  where  X  is  wavelength,  Lt(X) 
is  top-of-atmosphere  radiance,  L^X)  is  the  radiance  due  to  scattering  by  the  air  molecules  (Rayleigh  scattering) 

[10,1 1 ,1 2],  La(X)  is  the  radiance  due  to  scattering  by  aerosols,  L^X)  is  the  multiple  interaction  term  between  molecules 
and  aerosols,  t(X)  is  the  diffuse  transmittance  of  the  atmosphere  from  the  surface  to  the  sensor,  LwC(X)  is  the  radiance  due 
to  whitecaps  on  the  sea  surface  [13,14],  T(X)  is  the  direct  transmittance  from  the  surface  to  the  sensor,  Lg(X)  is  the 
specular  reflection  of  direct  sunlight  off  the  sea  surface  (sun-glitter)  [  1 5],  to(X)  is  the  diffuse  transmittance  of  the 
atmosphere  from  the  sun  to  the  surface,  0o  is  the  solar-zenith  angle,  and  [L*(X)]n  is  the  normalized  water-leaving 
radiance  (nl^)  due  to  photons  that  penetrate  the  sea  surface  and  are  backscattered  out  of  the  water  [16]. 

The  goal  of  atmospheric  correction  is  to  retrieve  nLw  accurately.  The  main  challenge  in  atmospheric  correction 
is  the  removal  of  Lpath(X)  from  Lt(X),  where  Lpath(X)  =  L^X)  +  La(X)  +  Lra(X).  Lpath(X)  contributes  about  90%  of  the  TOA 
radiance.  Lr(X)  can  be  removed  from  Lpjth(X)  by  using  standard  radiative  transfer  methods.  The  remaining  part  of 
Lpath(X),  La(X)  +  Lra(X),  is  estimated  from  Lt(X)  in  the  NIR  wavelengths.  We  assume  that  nLw  is  0  in  clear,  open-ocean 
regions  in  the  NIR  wavelengths.  However,  this  is  not  true  in  turbid,  coastal  areas.  By  assuming  that  nLw  is  0  in  the  NIR, 
this  provides  La(X)  +  Lra(X).  Based  on  La(X)  +  L„(X)  in  the  NIR,  an  estimate  is  made  of  La(X)  +  Lra(X)  in  the  visible. 

There  have  been  studies  conducted  over  the  past  decade  of  the  optical  properties  of  aerosol  types  [17,18,19,20,21,22,23]. 
This  has  lead  to  using  aerosol  modeling  to  extrapolate  from  the  NIR  to  the  visible  bands.  Figure  1  shows  these 
components  over  multiple  wavelengths  for  Gulf  of  Mexico,  May  4,  2010. 


Figure  1 .  MODlS-retrieved  Lt,  Ln  La,  and  nL^  for  412,  443,  488,  531,  547,  667,  748,  and  869  (nm)  wavelengths  at  the  AFRONET- 
OC  location  for  Gulf  of  Mexico,  May  4,  2010 

There  are  80  aerosol  models  (indexed  0  to  79)  to  choose  from  to  compute  La  4  Lra  during  atmospheric 
correction.  These  aerosol  models  are  used  for  computing  La(X)  4  Lra(X).  The  most  significant  digit  in  the  aerosol  model 
number  (ex:  the  6  in  model  64)  denotes  relative  humidity  index.  The  least  significant  digit  denotes  a  particle  size 
fraction  index.  Table  1  lists  the  different  relative  humidity  and  size  fraction  percentages,  along  with  their  corresponding 
aerosol  model  index.  Regarding  the  size  fractions,  a  size  fraction  of  20%  denotes  20%  fine  mode  and  80%  coarse  mode. 
For  example,  model  64  corresponds  to  a  relative  humidity  of  90%  and  a  size  fraction  of  20%,  with  20%  being  fine  mode 
and  80%  being  coarse  mode. 
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Table  1 .  Relative  humidity  and  size  fraction  percentages  and  corresponding  model  index,  models  ranging  from  0  to  79 


Relative  Humidity 

Index  (Most  Significant 
Digit  in  Model  Index) 

Relative  Humidity 
Percentage 

Size  Fraction  Index 
(Least  Significant  Digit 
in  Model  Index) 

Size  Fraction 

Percentage 

0 

30% 

0 

95% 

1 

50% 

1 

80% 

2 

70% 

2 

50% 

3 

75% 

3 

30% 

4 

80% 

4 

20% 

5 

85% 

5  • 

10% 

6 

90% 

6 

5% 

7 

95% 

7 

2% 

8 

1% 

9 

0% 

During  standard  atmospheric  processing,  there  are  two  aerosol  models  chosen  to  bound  La(A.)  +  L^A.)  in  the 
NIR.  The  aerosol  models  are  chosen  by  determining  which  two  aerosol  models  bound  e(748,869)  the  tightest,  where 
e(748,869)  is  a  ratio  for  wavelengths  748/869  [24].  This  ratio  is  used  to  select  bounding  aerosol  models.  The  relative 
humidity  index  is  calculated  before  e(748,869)  is  calculated,  so  e(748,869)  is  used  to  select  the  size  fraction  index. 
Once  the  two  bounding  aerosol  models  are  chosen  based  on  e(748,869) ,  interpolation  is  performed  between  the  two 
models  and  then  La(X)  +  Lrd(X)  is  retrieved  from  the  aerosol  model  lookup  tables. 


To  examine  the  effect  of  e(748,869)  and  model  selection  on  retrieved  nLw,  we  tested  all  80  aerosol  models 
individually  for  a  single  clear  day,  January  9,  2009,  in  Venice.  Instead  of  using  two  bounding  aerosol  models,  we  used  a 
fixed  aerosol  model  to  study  the  effects  of  all  80  aerosol  models  on  retrieved  nLw.  The  results  are  displayed  in  Figure  2. 
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Figure  2.  MOD1S  retrieved  epsilon  and  nL*(412)  for  all  80  aerosol  models  at  the  AERONET-OC  location  for  Venice,  Jan  9,  2009 


Figure  2  shows  that  relative  humidity  has  a  low  impact  on  retrieved  nLw  values.  To  describe  the  impact  on 
relative  humidity  and  size  fraction  on  retrieved  nLw  values,  consider  an  example  using  the  values  in  Figure  2  where  the 
bounding  aerosol  models  selected  are  14  (relative  humidity  of  50%,  size  fraction  of  20%)  and  15  (relative  humidity  of 
50%,  fine-mode  fraction  of  10%).  Model  14  and  15  depict  nLw(4I2)  ~  0.45  and  nLw(4I2)  ~  0.81,  respectively.  There  is 
significant  interpolation  between  these  two  models  to  derive  the  actual  nLw(412)  value.  This  is  due  to  the  distance  in 
nLw(4I2)  values  because  of  the  size  fraction,  not  the  relative  humidity.  If  the  model  selection  is  off  by  I  or  2  model 
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indices  with  regard  to  the  size  fraction,  then  there  are  large  errors  associated  with  the  retrieval  of  ,,1^(4 12).  Even  if  the 
models  are  correctly  selected,  the  interpolation  between  the  two  models  can  still  lead  to  large  errors  in  ,,1^(4 12). 
However,  the  relative  humidity  can  be  incorrectly  selected,  and  as  long  as  the  size  fraction  is  correct,  the  retrieved 
nLw(412)  value  should  still  match  closely  to  the  in  situ  value. 

To  further  demonstrate  how  size  fraction  is  a  larger  factor  than  relative  humidity'  in  retrieved  nLw  values, 
suppose  we  now  have  bounding  aerosol  models  74  and  75  from  Figure  2  instead  of  models  14  and  15  (relative  humidity 
changes  from  50%  to  95%  and  size  fraction  stays  the  same).  For  models  74  and  75,  ,,1^(412)  ~  0.4  and  nLw(41 2)  ~  0.6, 
respectively.  In  the  example  in  the  preceding  paragraph,  bounding  models  14  and  1 5  produce  nLw(412)  -  0.63  (halfway 
between  models  14  and  15).  Halfway  between  models  74  and  75,  nLw(412)  ~  0.5.  The  relative  humidity'  could  be  off  by 
several  model  indices  (in  this  case  6  models,  in  regards  to  the  relative  humidity  index),  and  the  nLw(41 2)  error  ~  0. 1 3.  To 
illustrate  the  impact  of  the  size  fraction,  now  suppose  the  bounding  aerosol  models  are  15  and  16  instead  of  14  and  15. 
For  models  15  and  16,  nLw(412)  ~  0.81  and  nL^(412)  ~  1.09,  respectively.  nLw(412)  interpolated  midway  between  these 
two  models  produces  nLw(412)  ~  0.95.  That  is  a  difference  of  0.32  from  models  14  and  1 5.  This  shows  that  if  incorrect 
aerosol  models  are  chosen  in  relation  to  the  size  fraction,  there  can  be  significant  errors  in  the  retrieved  nL^  values,  but 
smaller  errors  will  result  from  incorrect  model  selection  based  on  relative  humidity. 

3.  AUTOMATIC  VS  OPTIMAL  AEROSOL  MODEL  SELECTION 

Using  the  aerosol  models  selected  automatically  from  the  SeaWiFS/MODlS  atmospheric  correction  algorithm 
[25],  we  derive  nLw  from  MOD1S  1km  resolution  imagery  at  the  locations  of  the  AERONET-OC  sites  (Martha's 
Vineyard,  Venice,  and  Gulf  of  Mexico).  We  compare  these  satellite  values  to  the  AERONET-OC  measurements 
(matchups  within  3  hours,  at  least  50%  valid  pixels  in  a  5  x  5  box  centered  around  the  in  situ  pixel).  We  then  reprocess 
the  satellite  imagery  using  all  80  aerosol  models.  Rather  than  standard  processing  automatically  selecting  two  bounding 
aerosol  models,  a  single  aerosol  model  is  used  when  determining  the  optimal  model.  The  processing  is  repeated  80 
different  times,  once  for  each  model.  After  reprocessing  with  all  80  individual  aerosol  models,  we  again  compare  the 
derived  nLv,  products  to  the  in  situ  measurements.  The  model  yielding  the  closest  nU,  value  to  the  in  AERONET-OC 
value  is  selected  as  the  “optimal”  model.  For  each  scene  in  the  data  set,  we  are  only  interested  in  the  5  x  5  box  centered 
around  the  in  situ  pixel. 

For  the  analysis  of  Martha's  Viney  ard,  2010,  we  determine  the  optimal  aerosol  model  for  nLw  at  the  AERONET- 
OC  location  at  selected  visible  wavelengths.  Figure  3  displays  the  nLw  matchups  for  satellite  standard  processing  vs  in 
situ,  as  well  as  satellite  processing  with  optimal  aerosol  model  selection  vs  in  situ.  The  errors  (calculated  as  the  absolute 
value  of  ((MOD1S  -  AERONET-OC)  /  AERONET-OC)*  100)  are  greatly  reduced  with  the  optimal  model  selection. 
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Figure  3.  Automatic  aerosol  model  selection  (left  image)  vs  optimal  aerosol  model  selection  (right  image)  used  to  produce 
nLw  at  41 2,  443, 488,  547,  and  667  wavelengths  for  Martha’s  Vineyard,  2010 
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Using  the  same  data  set  seen  in  Figure  3,  in  Figure  4  we  show  the  optimal  aerosol  model  selected  for  computing 
nLw(41 2),  as  well  as  the  bounding  aerosol  models  selected  during  standard  atmospheric  processing.  Figure  4  shows  a 
wide  spread  of  optimal  aerosol  model  selections  throughout  2010  for  Martha's  Vineyard.  This  is  due  to  relative  humidity 
index  vary  ing  throughout  the  optimal  models  chosen  for  this  data  set.  The  optimal  model  for  a  particular  sample  in  the 
data  set  is  based  on  how  close  nLv,  is  to  the  AERONET-OC  value.  Because  of  this,  relative  humidity  is  not  really  used 
since  size  fraction  plays  a  more  important  role  in  retrieved  nL*  values.  This  does  not  mean  that  the  optimal  model  is 
stating  that  the  relative  humidity  is  incorrect  during  standard  atmospheric  correction.  It  also  does  not  mean  that  there  are 
not  other  models  that  are  capable  of  producing  values  that  closely  match  the  AERONET-OC  values.  For  example,  in 
Fig.  2  we  can  see  that  models  36,  46,  56,  and  66  all  yield  similar  nLw  values  (~  0.9). 


During  standard  atmospheric  correction,  the  modmin  bounding  aerosol  model’s  index  is  supposed  to  be  one 
higher  than  the  modmax  bounding  aerosol  model.  This  is  not  always  the  case.  There  are  four  points  in  Figure  4  where 
modmin  and  modmax  are  the  same.  They  are  the  5  x  5  box  means  centered  around  the  AERONET-OC  station  for  days 
80,  162,  23 1 ,  and  240.  This  is  the  case  when  e(748,869)  does  not  fall  between  the  lookup  tables  for  two  bounding 
aerosol  models.  When  modmin  and  modmax  are  the  same,  it  means  that  e(748,869)  is  either  lower  than  the  e(748,869) 
in  the  lookup  tables  for  the  lowest  available  bounding  aerosol  model  or  higher  than  the  €(748,869)  in  the  lookup  tables 
for  the  highest  available  bounding  aerosol  model.  When  this  occurs,  it  almost  always  leads  to  a  poor  estimate  of  the 
aerosol  composition,  which  in  turn  results  in  an  erroneous  estimate  of  nLw  values. 
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Figure  4.  Aulomalic  aerosol  selection  vs  oplimal  aerosol  model  seleclion  for  nLw(412)  for  Martha's  Vineyard,  2010 


Figure  5  uses  the  same  data  set  used  in  Figures  3  and  4.  In  Figure  5,  we  demonstrate  how  well  automatic 
aerosol  selection  (modmin  and  modmax  bounding  aerosol  models)  produces  nLw(412)  values  when  the  size  fraction 
matches  the  size  fraction  of  the  optimal  aerosol  model.  We  also  show  results  for  when  the  size  fraction  for  the  bounding 
aerosol  models  do  not  match  the  optimal  model.  The  nLw(41 2)  values  are  matched  against  AERONET-OC  nLw(4 1 2) 
values.  For  each  data  point,  we  are  comparing  how  well  the  size  fraction  (we  are  only  looking  at  size  fraction  here  and 
not  relative  humidity)  is  automatically  selected  compared  to  the  optimal  size  fraction  for  each  data  point  in  Figure  4. 
Because  the  automatic  aerosol  selection  uses  two  bounding  aerosol  models,  and  we  only  found  a  single  optimal  aerosol 
model  for  each  data  point,  we  must  determine  how  close  the  optimal  size  fraction  is  to  the  size  fractions  computed  for 
the  bounding  aerosol  models  during  the  automatic  aerosol  selection  process  of  standard  atmospheric  correction.  In 
Figure  5,  all  of  the  data  points  used  for  Martha’s  Vineyard  2010  are  separated  into  four  groups.  The  first  group  (modmax 
<=  optimal  <=  modmin)  is  for  any  data  points  where  the  optimal  size  fraction  is  in  between  the  bounding  aerosol 
models.  This  means  that  the  size  fraction  is  likely  to  have  been  correctly  selected  during  standard  atmospheric 
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correction.  For  example,  day  9  has  an  optimal  size  fraction  of  2,  a  5  x  5  box  average  for  modmin  of  2.96,  and  a  5  x  5 
box  average  for  modmax  of  1 .96.  To  clarify,  individual  points  in  the  scene  do  not  have  a  fractional  aerosol  model  index. 
For  modmax  to  have  a  value  of  1 .96,  24  of  25  points  in  the  5  x  5  box  have  a  value  of  2,  and  1  of  25  points  in  the  5  x  5 
box  have  a  value  of  1 .  The  second  group  (modmax  -  1  <=  optimal  <=  modmin  +  I )  is  for  data  points  where  the  size 
fractions  for  the  bounding  aerosol  models  do  not  fit  in  the  first  group,  but  are  only  plus  or  minus  I  size  fraction  away 
from  the  optimal  size  fraction.  For  example,  day  1 0  has  an  optimal  size  fraction  of  3,  a  5  x  5  box  average  for  modmin  of 
4.44,  and  a  5  x  5  box  average  for  modmax  of  3.44.  Day  1 0  belongs  in  the  second  group  because  (3.44  -  I  <=  3  <=  4.44 
+  1 ).  The  third  group  (modmax  -  2  <=  optimal  <=  modmin  +  2)  contains  points  +-  3  size  fraction  indices,  excluding 
those  in  groups  one  and  two.  The  fourth  group  (modmax  -  3+  <-  optimal  <=  modmin  +  3+)  is  for  all  data  points  that  do 
not  fit  in  groups  one,  two,  or  three. 

For  Figure  5,  we  remove  the  data  points  in  the  first  group  where  the  bounding  aerosol  models  have  the  same 
size  fraction  (modmin  and  modmax  are  both  equal  to  zero).  These  are  treated  as  bad  data  points,  for  reasons  previously 
discussed.  There  are  1 1  data  points  in  the  first  group,  1 2  in  the  second  group,  6  in  the  third  group,  and  7  in  the  fourth 
group.  This  means  that  standard  processing  accurately  selected  bounding  aerosol  models  for  only  1 1  of  the  total  36 
points.  The  average  error  for  the  data  points  in  the  first  group  is  1 9.26%,  37.85%  for  the  second  group,  35.60%  for  the 
third  group,  and  50.4%  for  the  fourth  group.  This  shows  that  if  the  bounding  aerosol  models’  fine  fraction  is  chosen 
correctly,  the  nLw  estimates  (nLw(412)  in  this  case)  are  reasonably  accurate.  If  the  size  fraction  is  off  by  just  one,  in  either 
direction,  it  can  greatly  impact  retrieved  nLw  values. 
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Figure  5.  Size  fraction  sensitivity  analysis  for  automatic  aerosol  selection  (size  fraction)  vs  AERONET-OC  fornLw(412)  for 
Martha's  Vineyard,  2010 


To  examine  model  variability,  in  Figure  6  we  show  the  ensemble  of  retrived  nLw  values  from  all  80  models  from 
the  MODIS  image  on  25  July  2010  at  the  Martha’s  Vineyard  AERONET-OC  site,  compared  to  the  in  situ  measurement. 
During  standard  automated  processing,  the  5  x  5  box  mean  surrounding  the  AERONET-OC  site  is  1 .07344,  and  the  two 
bounding  aerosol  models  are  54  and  55.  The  optimal  model  is  model  70,  which  produces  an  nLw(41 2)  of  0.7298,  which 
is  almost  identical  to  the  in  situ  value  of  0.7254.  Standard  processing  shows  a  relative  humidity  of  85%  (since  the  model 
index  is  in  the  50s),  and  the  optimal  model  shows  a  relative  humidity  of  95%  (since  the  model  index  is  in  the  70s).  The 
significant  difference  in  the  bounding  models  selected  during  standard  processing  and  the  optimal  model  is  the  size 
fraction.  The  modmax  bounding  model  (model  54)  produces  a  value  of  0.9222.  The  modmin  bounding  model  (model 
55)  produces  a  value  of  1 . 1 973.  We  can  take  the  mean  between  these  two  models  to  produce  the  standard  processing 
nLw(4l  2)  value  of  1 .07344.  Despite  the  optimal  aerosol  model  being  70  in  this  example,  there  are  multiple  aerosol 
models  that  can  select  an  nLw  value  that  matches  closely  to  the  AERONET-OC’s  measurement.  If  we  stay  in  the  standard 
processing  relative  humidity  (85%),  models  52  and  53  could  be  used  to  produce  a  better  nLw  value.  In  this  example, 
c(748,869)  is  computed  and  incorrect  bounding  aerosol  models  are  chosen.  Rather  than  models  54  and  55  being  chosen 
during  standard  atmospheric  correction,  if  models  52  and  53  were  chosen,  the  nLw  value  would  match  closely  to  the 
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in  situ  value.  This  is  an  example  of  how  during  standard  atmospheric  correction,  whenever  MOD1S  values  are  10%  ofT 
or  more  than  the  AERONET-OC  values,  it  is  usually  because  the  bounding  aerosol  models  were  incorrectly  selected. 


■  nLw(412)  Radiance 
•  optimal  Radiance 


Figure  6.  MOD1S  retrieved  nLvX412)  for  all  80  aerosol  models  at  the  AERONET-OC  location  for  Martha’s  Vineyard,  July  25,  2010 

For  the  analy  sis  of  Venice,  20 1 0,  we  determine  the  optimal  aerosol  model  at  nLw(4 1 2)  for  53  individual  points 
(each  point  is  a  separate  MOD1S  satellite  image;  each  individual  image  represents  a  single  day)  at  the  AERONET-OC 
location.  We  then  use  that  optimal  model  for  the  remaining  w  avelengths,  instead  of  calculating  the  optimal  model  for 
each  wavelength  as  we  did  with  the  Martha's  Vineyard  evaluation.  Figure  7  is  a  composite  figure  of  the  nL*  results.  The 
images  in  Figure  7  show  that  using  the  optimal  model  for  nLw(412)  has  good  performance  in  the  remaining  visible 
wavelengths. 


MODIS  nLw  Radiance  <mW/cm2/pm/sr)  MOOIS  nlw  Radiance  (mW/cm^m/sr) 

(5x5  Box  Mean)  (5x5  Box  Mean) 


Figure  7.  Automatic  aerosol  model  selection  (left  image)  vs  optimal  aerosol  model  selection  (right  image)  used  to  produce 
nLw  at  4 1 2,  443,  488,  547,  and  667  wavelengths  for  Venice,  20 1 0 
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Figure  8  shows  the  results  from  a  similar  analysis  for  the  Gulf  of  Mexico  AERONET-OC  site.  However,  the 
Gulf  of  Mexico  does  not  have  as  many  valid  points  for  2010  as  Martha's  Vineyard  or  Venice.  This  is  due  to  a  large 
amount  of  cloud  coverage  during  the  year,  as  well  as  the  AERONET-OC  station  being  unavailable  for  a  few  months  of 
the  year  due  to  instrument  calibration.  One  area  of  interest  in  this  analysis  is  that  for  two  of  the  MODIS  images,  there 
are  no  good  aerosol  models  available  that  are  capable  of  producing  a  matchup  close  to  the  corresponding  in  situ  value. 
There  are  five  data  points  in  Figure  8  that  correspond  to  these  two  MODIS  images  (two  nL^(412)  values,  two  nLw(443) 
values,  and  one  nL^(547)  value).  Three  of  these  five  values  (one  nLw(412)  value,  one  nLw(443)  value,  and  one  nLw(488) 
value)  correspond  to  day  176.  For  day  176,  both  modmin  and  modmax  are  equal  to  forty,  meaning  these  bounding 
aerosol  models  selected  during  standard  atmospheric  correction  are  incorrectly  chosen.  In  this  case,  there  is  no  possible 
aerosol  model  that  can  be  used  to  produce  a  good  result  for  ,,1^(412)  for  day  176.  This  is  an  instance  where  the  MODIS 
image  has  sporadic  cloud  coverage,  as  well  as  haze.  During  processing,  not  enough  pixels  were  flagged,  thus 
considering  this  day  to  be  a  valid  day.  When  there  are  no  possible  aerosol  models  capable  of  producing  values  that 
closely  match  AERONET-OC  values,  it  is  likely  because  pixels  were  not  flagged  as  invalid  pixels  when  they  should  be. 
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Figure  8.  Automatic  aerosol  model  selection  (left  image)  vs  optimal  aerosol  model  selection  (right  image)  used  to  produce 
nLw  at  412, 443, 488,  547,  and  667  wavelengths  for  Gulf  of  Mexico,  2010 


4.  CONCLUSIONS 

In  this  paper,  we  have  compared  satellite  nLw  retrievals  to  in  situ  for  Martha’s  Vineyard,  Venice,  and  the  Gulf  of 
Mexico.  During  standard  automated  atmospheric  correction,  selection  of  inappropriate  or  erroneous  bounding  aerosol 
models  can  dominate  the  errors  in  the  satellite  estimation  of  nLw.  Errors  in  nLw  will  in  turn  affect  downstream  bio-optical 
properties.  When  bounding  aerosol  models  are  incorrectly  chosen  during  atmospheric  correction,  it  is  the  size  fraction, 
rather  than  the  relative  humidity,  that  has  the  most  impact  on  retrieved  nLw  values.  If  bounding  aerosol  models  are 
incorrectly  selected,  there  is  usually  another  set  of  bounding  aerosol  models  within  the  same  relative  humidity  index  that 
are  capable  of  correctly  estimating  aerosol  values. 


5.  FUTURE  WORK 

There  usually  exists  a  temporal  or  regional  trend  of  over-estimation  or  under-estimation  of  nLw  satellite  values. 
Where  a  trend  exists,  we  can  devise  a  method  for  developing  ensembles  based  on  a  sampling  of  bounding  aerosol 
models.  These  ensembles  will  be  used  in  creating  new  nLw  values  that  better  resemble  values  from  nearby  AERONET- 
OC  stations. 

We  also  need  to  expand  our  study  from  a  concentrated  5x5  box  mean  around  the  AERONET-OC  point  to  the 
entire  scene.  If  we  change  the  bounding  aerosol  models  used  during  atmospheric  correction,  we  would  like  to  determine 
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if  the  change  in  bounding  aerosol  models  can  be  used  across  the  entire  scene.  This  can  be  evaluated  by  comparing 
chlorophyll  examinations  at  the  AERONET-OC  site  to  other  chlorophyll  samples  taken  somewhere  else  in  the  scene, 
perhaps  during  a  cruise.  We  need  to  further  examine  wavelengths  other  than  412  since  the  chlorophyll  algorithms  use 
other  visible  wavelengths. 

Another  method  for  improving  nLw  estimates  is  to  expand  the  set  of  aerosol  models.  More  models  are  needed  in 
the  middle  of  the  size  fraction  range  to  reduce  interpolation  between  bounding  aerosol  models. 
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